Driven interfaces in disordered media: determination of universality classes from 

experimental data 



00 
OS 

o> 

T— I 
Oh- 

(N 
(N 



43 

o 



I 

c 

o 
o 



> 

o 

m 
o> 
o 

00 

o> 
C3 



-a 
c 

o 
o 



X 
S3 



Reka Albert, Albert-Laszlo Barabasi* 
Department of Physics, University of Notre-Dame, Notre-Dame, IN 46556 

Nathan Carle, Andrew Dougherty 
Department of Physics, Lafayette College, Easton, PA 18042 
(February 1, 2008) 

While there have been important theoretical advances in understanding the universality classes 
of interfaces moving in porous media, the developed tools cannot be directly applied to experiments. 
Here we introduce a method that can distinguish the isotropic and directed percolation universality 
classes from snapshots of the interface profile. We test the method on discrete models whose uni- 
versality class is well known, and use it to identify the universality class of interfaces obtained in 
experiments on fluid flow in porous media. 



Interface motion in a random medium presents us with 
an archetype problem, with direct impact on various phe- 
nomena in condensed matter physics, including fluid flow 
in porous media |Q, domain growth in disordered mag- 
nets H and flux lines in disordered superconductors ||. 
In particular, much attention has been focused on un- 
derstanding the morphological evolution of an interface 
driven through a disordered medium. Two-phase fluid 
flow experiments have provided evidence that the mor- 
phology of the interface can be either self-similar or self- 
affine Q] . The self-similar morphology has been success- 
fully described by various percolation based models, such 
as invasion percolation The motion and roughening 
of self-afnne interface morphologies can be quantified by 
the global interface width w(L,t), where L is the system 
size. The study of discrete models and continuum growth 
equations lead to the observation that the width follows 



w(L,t) 



t 13 if t < t x 
L a if t > t x 



(1) 



where (3 is the growth exponent and a is the roughness 
exponent. While both experiments and models have con- 
firmed the applicability of scaling concepts, the obtained 
roughness exponents are scattered between 0.6 and 1.25, 
H questioning the existence of universality, the founda- 
tion of the scaling hypothesis (0) . 

Motivated by the success of the Kardar-Parisi-Zhang 
(KPZ) equation in describing interface motion with ther- 
mal noise, it has been proposed that interfaces in porous 
media are described by the quenched KPZ (QKPZ) equa- 
tion 

^L=F + vV 2 h+~(Vh) 2 +r,(x,h), (2) 

where rj(x., h) represents the quenched noise in the 
medium. This equation predicts the existence of a depin- 
ning transition: for driving forces F < F c the interface is 
pinned, and for F > F c it moves with a velocity v ~ f e , 
where 9 is the velocity exponent and f = (F — F c )/F c . 



Numerical R| and analytical JsJ studies indicate that 
the discrete models can be grouped into two universality 
classes depending on the behavior of the QKPZ nonlin- 
earity A. Isotropic models (i.e. models that have no 
growth direction determined by the random medium) 
have A = or A - > as / — > 0. The scaling expo- 
nents of the isotropic universality class can be deter- 
mined analytically from (||) with A = 0, obtaining on = 1, 
Pi — 0.75 and Oi = 0.33 in one dimension ||. However, 
anisotropy in the medium can induce a A term that di- 
verges at the depinning transition. The exponents char- 
acterizing this anisotropic universality class can be ob- 
tained from an exact mapping to directed percolation, 
providing a a = 0.63, (3 a = 0.63 and 6 a = 0.63 0,0. 

While in numerical simulations the universality class 
can be identified by measuring A for a tilted interface 
0, this method cannot be applied experimentally, since 
a constant tilt cannot be sustained. The scaling expo- 
nents obtained for paper wetting JlC],|l2| are in excellent 
agreement with the prediction of the anisotropic univer- 
sality class. However, experiments on fluid flow between 
glass beads provided exponents between 0.6 and 0.9 |jl3| , 
not allowing to identify the universality class to which 
they belong. Thus, in order to properly characterize the 
experimental results we need to develop methods which, 
without relying on the determination of the scaling ex- 
ponents, can identify and distinguish the two universality 
classes. Here we make a major step in this direction by 
introducing a method that can identify the universality 
class from snapshots of the interface at different time in- 
tervals. After demonstrating that the method succesfully 
identifies the universality class of well known models, we 
apply it to identify the universality class of interfaces ob- 
tained in experimental investigations. 

Consider an interface h(i,t) moving in a porous 
medium. Since the interface is rough, we can define a 
local slope Sj as being the linear fit to the interface in the 
region (i — 8/2, i + 5/2) (see Fig. 1). If there are slope 
dependent terms governing the motion of the interface, 



(such as the \(Vh) term in (|2j)), the local velocity Ui is 
expected to depend on the local tilt Si. 
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The main hypothesis of this paper is that the local 
slope- dependent velocities satisfy similar scaling laws to 
the average velocity. We expect that the local velocity 
u(F, s), corresponding to the local slope s, obeys the scal- 
ing law 



u(F,s)~(F-F^s)) t 



(3) 



where u(F, s) is the average local velocity of interface 
segments with slope s, i^(s) is the depinning threshold 
corresponding to these segments, and 9' is the local ve- 
locity exponent. 

In an isotropic medium there is no selected growth di- 
rection, thus the depinning threshold and the velocity 
exponent are independent of slope, i.e. F' c (s) = F c and 
9' = 9. Consequently, the local velocities have the form 
u(F, s) = w(s)(F — F c ) e , analogous to the scaling of the 
average velocity v(F) — v (F — F c ) e , i.e. the only slope 
dependence comes in the non-universal prefactor w(s). 
Thus we expect that for an isotropic medium the ratio 



Vi(F,s) 



u(F, s) w(s) 



(4) 



is independent of the driving force F. 

In contrast, one of the distinguishing features of an 
anisotropic medium is that the depinning threshold, 
.Fj!(s), decreases with s. Furthermore, it has been shown 
H that for globally tilted interfaces (i.e. interfaces which 
have a nonzero global slope m), the velocity exponent 
is different from 9 a — 0.63 obtained for an untilted in- 
terface, and it is equal to 1. Consequently, we cannot 
exclude the possibility that for anisotropic media 9' a de- 
fined in (g) is different from 9 a . Thus for the anisotropic 
universality class the ratio V a = u(F, s)/v(F) has a non- 
trivial F -dependence: 



Va(F,s) 



w(s)(F - 

v (F-F c )^ 



(5) 



In conclusion, we expect that if V{F, s) is independent 
of the driving force, the interface belongs to the isotropic 
universality class, while a systematic F dependence is an 
indication that the interface belongs to the anisotropic 
universality class. Note that in experiments the driving 
force is not available, but we can replace it with the av- 
erage velocity, because v(F) ~ (F — F c ) e . 

The advantage of this method is that it can be applied 
directly to both experimental and numerical data, as long 
as we have snapshots of the interface at not too distant 
time intervals. The method is applied as follows: a dis- 
cretized interface of total length L is partitioned into seg- 
ments of length 6. The local slope s(i) for each of the L/6 
segments is determined by fitting a line to it in the region 
(i-6/2,i + 6/2), where i = 6/2, 36/2, .., L-6/2. We then 
repeat the same partitioning for the interface captured at 
time t + r, and calculate the local velocities of each of the 
segments centered at i as u(i, s) = [h(i,t + r) — h(i, t)]/r. 
As a first step we calculate the average velocity of all 



segments with slope s when the average velocity of the 
interface is v, u{v,s) = (1/N(s)) $^ti(i, s), where the 
sum goes over all segments which have the same slope 
s, and N(s) is the number of such segments (naturally, 
^ s N(s) = L/6). If u(v, s) plotted as a function of s is 
a parabola, it indicates that a slope dependent nonlinear 
term is present in the growth equation. Finally, plotting 
V(v, s) = u(v, s)/v allows us to distinguish the two uni- 
versality classes: if the u(v, s) curves collapse into a single 
one, we expect the model to belong to the isotropic uni- 
versality class, while if a systematic v dependence persists 
after rescaling, it belongs to the anisotropic universality 
class. 

To test the applicability of the method, we calculate 
V(v,s) for models whose universality class is known ||, 
namely the Random Field Ising Model (RFIM) belong- 
ing to the isotropic universality class pj], and the di- 
rected percolation depinning (DPD2) model |p|,pl|, that 
is known to belong to the anisotropic universality class. 
These simulations also offer support for the scaling hy- 
pothesis (||). 

Random Field Ising model — The definition and the 
simulation techniques of the RFIM are described in de- 
tail in Ref. We have found that the depinning 
threshold of the local segments is independent of their 
slope and it is equal to the depinning threshold of the 
average velocity, F c . Also, we found that the velocity ex- 
ponent of the local velocities is constant and it is equal to 
6>;. Figure 2a shows the local velocities u(v, s) versus the 
local slope s for different driving forces. The parabolic 
shape of the curves indicates that the KPZ nonlinearity 
is present, in agreement with the results obtained from 
the global tilt measurements 0. As Fig. 2b indicates, 
the V(v, s) curves collapse into a single function, thus 
confirming that the RFIM belongs to the isotropic uni- 
versality class. 

DPD model— The DPD2 model, introduced by Tang 
and Leschhorn, is described in detail in Ref. 0,0- As 
a first step, we investigated the slope-dependence of F' c 
and 9' . The scaling of u(F, s) has the following features: 

(a) F' c (s) decreases with |s| obeying the same scaling 
law as F c (m) corresponding to tilted interfaces B; 

(b) The velocity exponent is independent of slope and 
it has the value 9' a — 1. Consequently, while the average 
velocity has an exponent 9 a = 0.63, the local segments 
are all tilted from their local hard direction of depinning, 
and thus they have 9' a — 1 M. 

As shown in figure 3a the parabolic shape of the u(v, s) 
curves indicates the presence of a A(V/i) 2 term. However, 
in contrast with the results obtained for the RFIM, there 
is a systematic shift in the V a (v, s) curves with increasing 
v, confirming that the model belongs to the anisotropic 
universality class. 

Experiments — The experiments were performed in a 
thin rectangular plexiglass cell of 125 x 41 x 0.16 cm, sim- 
ilar to that used in Ref. jig] , and filled with unconsoli- 
dated glass beads approximately 170±20 /im in diameter. 
The apparatus was continually agitated during filling to 
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help ensure a reasonably uniform packing of the beads. 
Water was injected by a syringe pump through a series 
of small holes in a tube aligned along one of the narrow 
ends of the cell, approximately mimicking the ideal case 
of uniform injection along a line. Images were recorded at 
2-second intervals with a CCD camera and digitized with 
a resolution of 512 x 480 pixels. Only the central 17 cm of 
the system was recorded to avoid edge effects. The overall 
scale was 270 /zm/pixel, or about 1.5 beads/pixel. All ve- 
locities were measured in units of bead diameters/second. 

Fig. 4a shows the local velocity curves obtained from 
the experimental data. Their shape is consistent with a 
par abola, indicating the presence of the KPZ nonlinearity 
jig . Fig. 4b presents V(F, s) versus the local slopes for 
six different experiments. While there is some scattering 
in the data collapse due to the asymmetry of the curves 
p6[ , in general we can see that there is no systematic shift 
with the increasing velocity, indicating that while there is 
a KPZ nonlinearity (suggested by the parabolic shape of 
V), the experiments belong to the isotropic universality 
class. 

A convincing summary of our results is shown in Fig. 5, 
where we plot V(v,0) as a function of v for the models 
corresponding to the two universality classes. We see that 
for the RFIM Vi(v, 0) is independent of the average veloc- 
ity, while for the DPD model, V a (v, 0) increases systemat- 
ically with v. Indeed, the driving force dependence of the 
s = portions is (F - F c ) e '- e . In the case of the DPD2 
model 9 a = 0.63 and Q' a = 1, thus V a {v,Q) increases with 
F - F c . On the other hand, V a (v,0) ~ v'- 1 '^^" in- 
creases with v. On the same plot we show V(v,0) de- 
termined from the experimental data, which, within the 
error bars, is consistent with the line corresponding to the 
RFIM and deviates considerably from the DPD results, 
underlying our conclusion that the experiments belong 
to the isotropic universality class. 
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FIG. 1. In the interval of length S around the position i 
the interface is fitted with a straight line of slope Si. Assum- 
ing that in a consecutive interface h(i, t + r) the slope is still 
Si, the local velocity of point i is Ah/r. 

FIG. 2. RFIM model, (a) Local velocity u(v, s) versus the 
local slope s for driving forces / = (H — H c )/ H c ranging from 
0.015 (bottom curve) to 0.076 (top curve). The system size 
is L — 500 and each result was averaged over 200 realizations 
of the disorder, (b) The same curves rescaled by the average 
velocities. 

FIG. 3. DPD2 model, (a) u(v, s) versus s for driving forces 
/ = (F - F c )/F c ranging from 0.016 (bottom curve) to 0.07 
(top curve). The system size is L — 20000 and each result 
was averaged over 500 realizations of the disorder, (b) The 
same curves rescaled by the average velocity. 
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FIG. 4. (a) Dependence of u(v, s) on s for the interfaces 
obtained in the experiments. The average velocity is given in 
bead diameters/second, (b) The same curves rescaled by the 
average velocities. 

FIG. 5. Comparison between the behavior of V(v, 0) with v 
for the two universality classes. Both variables were rescaled 
by v m i„. While Vi(v,0) is constant, V a {v,0) increases with 
v. On the same plot we show the values obtained from the 
experimental interfaces. The error bars for the experimental 
data were derived from the fluctuations of the average velocity 
of the interface. 
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